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Abstract 

A measurement of the CP asymmetry in — )■ K*^^'^fi~ decays is presented, based 
on 1.0 fb~^ of pp cohision data recorded by the LHCb experiment during 2011. The 
measurement is performed in six bins of invariant mass squared of the pair, 
excluding the J/ip and tp{2S) resonance regions. Production and detection asymme- 
tries are removed using the — t- J/i/jK*^ decay as a control mode. The integrated 
CP asymmetry is found to be —0.072 it 0.040 (stat.) ± 0.005 (syst.), consistent with 
the Standard Model. 
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The decay — )■ K*^{^ K^n )fi^iJ, is a flavour changing neutral current process which 
proceeds via electroweak loop and box diagrams in the Standard Model (SM) (l). The 
decay is highly suppressed in the SM and therefore physics beyond the SM such as su- 
persymmetry f2| can contribute with a comparable amplitude via gluino or chargino loop 
diagrams. A number of observables are sensitive to such contributions, including the 
partial rate of the decay, the fi~^^~ forward-backward asymmetry (Afb) and the CP 
asymmetry (Acp)- The CP asymmetry for — )■ K*^fi~^^~ is defined as 

where F is the decay rate and the initial flavour of the B meson is tagged by the charge 
of the kaon from the K* decay. The CP asymmetry is predicted to be of the order 10~^ 
in the SM [3|[4], but is sensitive to physics beyond the SM that changes the operator 
basis by modifying the mixture of the vector and axial- vector components (5[[6|. Some 
models that include new phenomena enhance the observed CP asymmetry up to ±0.15 |7|. 
The theoretical prediction within a given model has a small error as the form factor 
uncertainties, which are the dominant theoretical errors for the decay rate, cancel in the 
ratio. 

The CP asymmetry in B^ — )■ K*^fi^fi~ decays has previously been measured by the 
Belle [s] and BaBar [o] collaborations, with both results consistent with the SM. The 
LHCb collaboration has recently demonstrated its potential in this area with the most 
precise measurement of ApB flO], and in this Letter, the measurement of the CP asym- 
metry by LHCb is presented. 

The LHCb detector is a single-arm forward spectrometer covering the pseudo- 
rapidity range 2 < rj < 5, designed for the study of particles containing 6 or c quarks. 
The detector includes a high precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of 
silicon-strip detectors and straw drift-tubes placed downstream. The combined tracking 
system has a momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at 
100 GeV/c, and an impact parameter (IP) resolution of 20 [im for tracks with high trans- 
verse momentum (pt)- Charged hadrons are identified using two ring-imaging Cherenkov 
detectors. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and pre-shower detectors, an electromagnetic calorimeter 
and a hadronic calorimeter. Muons are identified by a system composed of alternating 
layers of iron and multiwire proportional chambers. The trigger consists of a hardware 
stage, based on information from the calorimeter and muon systems, followed by a soft- 
ware stage which makes use of a full event reconstruction. 

The simulated events used in this analysis are produced using the Pythia 6.4 genera- 



tor 12 , with a choice of parameters specifically configured for LHCb 13 . The EvtGen 



package 14 describes the decay of the particles and the Geant4 toolkit 15 simulates 



the detector response, implemented as described in Ref. 16 . QED radiative corrections 



are generated with the Photos package 17 
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The events used in the analysis are selected by a dedicated muon hardware trigger and 



then by one or more of a set of different of muon and topological software triggers 18,19 



The hardware trigger requires the muons to have pt greater than 1.48GeV/c, and the 
software trigger requires that one of the final state particles to have both > 0.8 GeV/c 

Triggered candidates 



10 . The first stage is a 



and IP with respect to all pp interaction vertices > 100 |^m 19 
are subject to the same two-stage selection as that used in Ref. 

cut-based selection, which includes requirements on the 5° candidate's vertex fit x^; the 
flight distance ^^id invariant mass and each track's impact parameters with respect to 
any interaction vertex, and polar angle. Background from misidentified kaon and pion 
tracks is removed using information from the particle identification (PID) system and 
muon tracks are required to have hits in the muon system. Finally, the production vertex 
of the candidate must lie within 5 mm of the beam axis in the transverse directions, 
and within 200 mm of the average interaction position in the beam (z) direction. 

In the second stage, the candidates must pass a multivariate selection that uses a 



boosted decision tree (BDT) [20] that implements the AdaBoost algorithm 21 . This is 



a tighter selection which takes into account other variables including the decay time and 
flight direction of the candidates, the of the hadrons, measures of the track and 
vertex quality, and PID information for the daughter tracks. For the rest of the Letter, 
the inclusion of charge conjugate modes is implied unless explicitly stated. 

In order to obtain a clean sample of B^ — )■ K*^ix^ ix~ decays, the cc resonant decays 
B^ — )■ J/ipK*^ and B^ — )■ tlj{2S)K*^ are removed by excluding events with fi~^fi~ invariant 
mass, m^+^-, satisfying 2.95 < m^+^- < 3.18GeV/c^ or 3.59 < m^+^- < 3.77GeV/c^. 
If rnK+TT-fi+fj.- < 5.23GeV/c^, then the vetoes are extended downwards by 0.15GeV/c^ to 
remove the radiative tails of the resonances. Backgrounds involving misidentified particles 
are vetoed using cuts on the masses of the B^ and K*^ mesons and the /i"*"//" pair, as 
well as using the PID information for the daughter particles. These include B^ 
candidates in which a kaon has been misidentified as a pion, 5° J/i)K*^ candidates 
where a hadron is swapped with a muon and B^ — )■ K^^^^~ candidates which combine 



with a random low momentum pion. The vetoes are described fully in Ref. 10 . A 
further background, which dilutes the level of Acp measured, is given by B^ — )■ K*^fi'^fi~ 
candidates in which the kaon and pion are misidentified as each other. However, the 
vetoes reduce the level of contamination to 0.8% of the B^ — )■ K*^n'^fi~ yield. All B^ 
candidates must have a mass in the range 5.15 — 5.80GeV/c^, the tight low mass edge of 
this window serves to remove background from partially reconstructed B meson decays. 
All K*^ candidates must have an invariant mass of the kaon-pion pair within 0.1 GeV/c^ of 
the nominal K*^{892) mass. A proton veto, using PID information from a neural network, 
is also applied to remove background from Af, decays, where a proton in the final state is 
misidentified as a kaon or pion in the B^ — )■ K*^fi'^fi~ decay. 

Approximately 2% of selected events contain two B^ — K*^^'^^~ candidates which 
have tracks in common. The majority of these candidates arise from swapping the as- 
signment of the kaon and pion hypothesis. As the charges of the kaon and pion tag the 
flavour of the B meson these duplicate candidates can bias the measured value of Acp- 
This is accounted for by randomly removing one of the two candidates from the sample. 
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This process is repeated many times over the full sample with a different random seed in 
each case and the average measured value of Acp is taken as the result. 

In order to obtain an accurate measurement of Acp the differences in the production 
rates (R) and detection efficiencies (e) between the 5° — t- K*^fi~^fi~ and — > K*^^^^~ 
modes must be accounted for. Assuming all asymmetries are small the raw measured 
asymmetry may be expressed as 

-4raw = -Acp + + Ab, (2) 



where the production asymmetry, which is of the order of 1% 22 , is de- 
fined as = [R (5°) - R (5°)] / [R (5°) + R (5°)] and the detection asymmetry is 
Ab = [e (-8°) — e (-8°)] / [e (5°) + e (-8°)] . The production asymmetry is diluted through 
B^-B^ oscillations by a factor n 

_ e(t)e~^* cos Amtdt 

where t, F, and Am are the decay time, mean decay rate, and mass difference between the 
light and heavy eigenstates of the B^ meson respectively. The quantity is dominated 
by the K^t:~ / K~tt^ detection asymmetry which arises due to left-right asymmetries in 
the LHCb detector and different interactions of positively and negatively charged tracks 
with the detector material. The left-right asymmetry is cancelled by taking an average 
with equal weights of the CP asymmetries measured in two independent data samples 
with opposite polarities of the LHCb dipole magnet. These data samples correspond to 
61% and 39% of the total data sample. 

The production and interaction asymmetries are corrected for using the B^ — )■ J /ipK*^ 
decay mode as a control channel. Since B^ — > K*^fi~^fi~ and 5^ — > J/ipK*^ decays have 
the same final state and similar kinematics, the measured raw asymmetry for B — t- J/ip K*^ 
decays may be simply expressed as ^raw — ^ J /ipK*^) = kAp + Ab, assuming the 
CP asymmetry in B^ — )■ J/ipK*^ decays is small. B^ — )■ J/ipK*^ proceeds via a 6 — )■ 
CCS transition, as does the decay mode B^ — J/ipK^, and hence should have a CP 
asymmetry similar to Acp{B+ J/ipK+) = (1 ± 7) x 10"^ |23|[24). The CP asymmetry 
in B^ — )■ K*^fi~^fi~ decays is then calculated as 

Acp = Akaw (5° -> K*V^/i") - ^raw {B' J/i^K*^) . (4) 

Non-cancelling asymmetries due to differences between the kinematics of B^ — )■ K*^fi^fi~ 
and B^ — )• J /ipK*^ decays are considered as systematic effects. 

The full data sample, containing approximately 900 B^ — )■ K*^fi'^^~ signal decays, is 
split into the six bins of invariant mass squared (g^) used by the LHCb, Belle, and 

CDF angular analyses [8 , 10 ,25 . An additional bin of 1 < < QGeV'^/c'^ is used, to be 



compared to the theoretical prediction in Ref. [4]. The values of Acp are measured using a 
simultaneous unbinned maximum-likelihood fit to the B^ — )■ K*^fi'^fi~ and B^ — >■ J/ipK*^ 
invariant mass distributions in the range 5.15 — 5.80 GeV/c^. The simultaneous fit in each 
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bin spans eight data samples, split between the initial particles and the decay 
modes — )■ K*^fi~^fi~ and B° — )■ J/ipK*^, and magnet polarity, where the 5° — )■ J j^K*^ 
sample is common to all bins. This fit returns two values of Acp-, one for each magnet 
polarity, and an average with equal weights is made to find the value of Acp in that 
bin. An integrated value of Acp over all is also calculated. 

The signal invariant mass distributions for the B^ — )■ K*^fi^fi^ and B^ — )■ J/ipK*^ 
decays are modelled using the sum of two Crystal Ball functions 26 with common peak 
and tail parameters but different widths. The values of the tail parameters are determined 
from fits to simulated events and fixed in the fit. Combinatorial background arising 
from the random misassociation of tracks to form a B^ candidate is modelled using an 
exponential function. The B^ — J/ipK*^ fit also accounts for a peaking B^ — )■ J/ipK*^ 
contribution, which has the same shape as the signal and an expected yield which is 
(0.7 ± 0.2)% of that of 5° J/ipK*^ |27]. In the simuhaneous fit, the signal shape is 
the same for the two modes, but the signal and background yields and the exponential 
background parameter may vary. Figure [l] shows the mass fit to the B^ — )■ K*^^'^^~ 
decay in the full range. 

Many sources of systematic uncertainty cancel in the difference of the raw asymmetries 
between B^ — )■ K*'^fi~^fi~ and B^ — )■ J/ipK*^ decays and in the average of CP asymme- 
tries measured using data recorded with opposite magnet polarities. However, systematic 
uncertainties can arise from residual non-cancelling asymmetries due to the different kine- 
matic behaviour of B'^ — )■ K*^fi'^fi~ and B^ — )■ J/ipK*^ decays. The effect is estimated 
by reweighting B^ — J/ tpK*^ candidates so that their kinematic variables are distributed 
in the same way as for 5° — )■ _ft'*°yU+/i" candidates. The value of Ajiaw — ?• J /ipK*^) 
is then calculated for these reweighted events and the difference from the default value is 
taken as the systematic uncertainty. This procedure is carried out separately for a num- 
ber of quantities including the p, px, and pseudorapidity of the 5° and the K*^ mesons. 
The total systematic uncertainty associated to the different kinematic behaviour of the 
two decays is calculated by adding each individual contribution in quadrature. This is 
conservative, as many of the variables are correlated. 

The random removal of multiple candidates discussed above also introduces a system- 
atic uncertainty on Acp- The uncertainty on the mean value of Acp from the ten different 
random removals is taken as the systematic uncertainty. 

The forward-backward asymmetry in B^ — > K*^fi'^fi~ decays 10 , which varies as 



a function of g^, causes the positive muons to have different average momenta, with 
the opposite effect for the charge conjugate mode. Different detection efficiencies for 
positive and negative muons will introduce an asymmetry that cannot be accounted 
for by the B^ — )■ J/ipK*^ decay, as this mode does not have a comparable forward- 
backward asymmetry. The selection efficiencies for positive and negative muons are eval- 
uated using muons from J/ip decay in data and the resulting asymmetry in the selected 
B^ — > K*^fi'^fi~ sample is calculated in each bin. 

A number of possible effects due to the choice of model for the mass fit are considered. 
The signal model is replaced with a sum of two Gaussian distributions and a possible 
difference in the mass resolution for B^ — )■ K*^fi~^fi~ and B^ — )■ J/ipK*^ decays is inves- 
tigated by allowing the width of the B^ — K*^fi'^fi~ signal peak to vary in a range of 
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Figure 1: Mass fits for — > K*^jj,~^j2~ decays used to extract tlie integrated CP asym- 
metry. Tlie curves displayed are tlie full mass fit (blue, solid), the signal peak (red, 
short-dashed), and the background (grey, long-dashed). The mass fits on the top row 
correspond to the (a) and (b) B^ decays for one magnet polarity, while the bottom 
row shows the mass fits for (c) B^ and (d) B^ for the reverse polarity. 

0.7—1.3 times that of the B^ — )■ J/ipK*^ model. As a further cross-check, Acp is calcu- 
lated using a weighted average of the measurements from the six bins and the result is 
found to be consistent with that obtained from the integrated dataset. 

The results of the full Acp fit are presented in Table |2] and Figure [2] The raw asym- 
metry in B^ —7- J/ipK*^ decays is measured as 

^RAw [B^ -> J/^K*^) = -0.0110 ± 0.0032 (stat.) ± 0.0006 (syst.). 

The CP asymmetry integrated over the full range is calculated and found to be 

Acp (5° ^ K*^fi+fi-) = -0.072 ± 0.040 (stat.) ± 0.005 (syst.). 

The result is consistent with previous measurements made by Belle [s], Acp (B — )■ 
K*l+l-) = -0.10 ± 0.10 (stat.) ± 0.01 (syst.), and BaBar g, Acp {B K*l+l-) = 
0.03 ± 0.13 (stat.) ± 0.01 (syst.). This is significantly the most precise measurement of 
Acp in B° K*^ii^ii~ decays to date. 
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Table 1: Systematic uncertainties on Acp, from residual kinematic asymmetries, muon 
asymmetry, choice of signal model, and the modelling of the mass resolution, for each 
bin. The total uncertainty is calculated by adding the individual uncertainties in 
quadrature. 





Sources of systematic uncertainties 










multiple 


residual 


yU^ detection 


signal 


mass 




region ( < 




cands. 


asymmetries 


asymmetry 


model 


resol. 


Total 


0.05 < 




< 2.00 


0.002 


0.007 


0.005 


0.005 


0.001 


0.010 


2.00 < 




< 4.30 


0.006 


0.007 


0.006 


0.007 


0.010 


0.016 


4.30 < 




< 8.68 


0.004 


0.003 


0.006 


0.004 


0.003 


0.010 


10.09 < 




< 12.86 


0.003 


0.007 


0.009 


0.001 


0.002 


0.011 


14.18 < 




< 16.00 


0.001 


0.006 


0.007 


0.001 


0.001 


0.009 


16.00 < 




< 20.00 


0.003 


0.005 


0.003 


0.003 


0.009 


0.012 


1.00 < 




< 6.00 


0.001 


0.006 


0.005 


0.002 


0.003 


0.009 


0.05 < 




< 20.00 


0.002 


0.002 


0.005 


0.001 


0.001 


0.005 



Table 2: Values of Acp for — > K*^iJL^iJi in the six bins used in the analysis. 









signal 




statistical 


systematic 


total 


region (GeVYc"') 


yield 


Acp 


uncertainty 


uncertainty 


uncertainty 


0.05 < 




< 2.00 


168.2±14.6 


-0.196 


0.094 


0.010 


0.095 


2.00 < 


q' 


< 4.30 


72.1±10.7 


-0.098 


0.153 


0.016 


0.154 


4.30 < 


q' 


< 8.68 


266.3±19.1 


-0.021 


0.073 


0.010 


0.075 


10.09 < 


q' 


< 12.86 


156.7±15.1 


-0.054 


0.097 


0.011 


0.098 


14.18 < 


q' 


< 16.00 


116.2±12.1 


-0.201 


0.104 


0.009 


0.104 


16.00 < 




< 20.00 


128.2±12.9 


0.089 


0.100 


0.012 


0.101 


1.00 < 




< 6.00 


194.0±17.0 


-0.058 


0.064 


0.009 


0.064 


0.05 < 


q' 


< 20.00 


904.1±35.3 


-0.072 


0.040 


0.005 


0.040 
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Figure 2: Fitted value of Acp in — > K*^fi^fj,~ decays in bins of the invariant 
mass squared (g^). The red vertical lines mark the charmonium vetoes. The points are 
plotted at the mean value of in each bin. The uncertainties on each Acp value are the 
statistical and systematic uncertainties added in quadrature. The dashed hne corresponds 
to the integrated value, and the grey band is the la uncertainty on this value. 
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